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ABSTRACT 

We present an updated phenomenological galaxy evolution model to fit the Spitzer 24, 70 and 
160 /um number counts as well as all the previous mid and far infrared observations. Only a minor 
change of the co-moving luminosity density distribution in the previous model (Lagache, Dole, Puget 
2003), combined with a slight modification of the starburst template spectra mainly between 12 and 
30 /mi, are required to fit all the data available. We show that the peak in the MIPS 24 /im counts is 
dominated by galaxies with redshift between 1 and 2, with a non negligible contribution from the z > 
2 galaxies (~ 30% at S=0.2 mJy). The very close agreement between the model and number counts 
at 15 and 24 /mi strikingly implies that (1) the PAHs (Policyclic Aromatic Hydrocarbons) features 
remain prominent in the redshift band 0.5 to 2.5 and (2) the IR energy output has to be dominated 
by ~3 10 11 Lq to ~3 10 12 L Q galaxies from redshift 0.5 to 2.5. Combining Spitzer with the Infrared 
Space Observatory (ISO) deep cosmological surveys gives for the first time an unbiased view of the 
infrared Universe from z=0 to z=2.5. 
Subject headings: infrared: galaxies -galaxies: evolution - galaxies: model 



1. INTRODUCTION 

A Cosmic far-Infrared Background (CIB) that would 
trace the peak of the star-formation and metal produc- 
tion in galaxy assembly has long been predicted. The 
first observational evidence for such a background was 
reported by Puget et al. (1996) (see Hauser & Dwek 
2001, for a review). This discovery, together with recent 
cosmological surveys in the infrared (IR) and submillimc- 
ter, has opened new perspectives on our understanding of 
galaxy formation and evolution. The surprisingly large 
amount of energy in the CIB shows that it is crucial to 
probe the galaxies contributing to it to understand when 
and how the bulk of stars formed in the Universe. 
To understand the sources contributing to the CIB and 
to interpret the deep source counts, we have developed 
a phenomenological model that constrains in a simple 
way the IR luminosity function evolution with redshift, 
and fits all the existing source counts consistent with the 
redshift distribution, the CIB intensity, and, for the first 
time, the CIB fluctuation observations from the mid-IR 
to the submm range (Lagache et al. 2003). Spitzer has 
provided deep new multiwavelength source counts from 
24 to 160 um (Papovich et al. 2004, Dole et al. 2004a). 
We use these new data to search for a new optimisation of 
the Lagache et al. (2003) model parameters. A remark- 
able result is that only a minor change of the co-moving 
luminosity density distribution combined with a slight 
modification of the starburst model spectra, mainly for 
the 12 < A <30 /mi range, are required to fit all the 
previous data together with the new constraints. 

2. THE MODEL 

The model is discussed extensively in Lagache et al. 
(2003). Our goal is to build the simplest model with 
the fewest free parameters and separate ingredients. We 
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fix the cosmology to fl\=0.7, fio=0.3 and h=0.65. We 
assume that IR galaxies are mostly powered by star for- 
mation and hence we use Spectral Energy Distributions 
(SEDs) typical of star- forming galaxies 3 . Although some 
of the galaxies will have AGN-dominated SEDs, they 
are a small enough fraction (< 10%: Alonso Herrero et 
al. 2004) that they do not affect the results significantly. 
We therefore construct 'normal' and starburst galaxy 
template SEDs: a single form of SED is associated with 
each activity type and luminosity. We assume that the 
Luminosity Function (LF) is represented by these two 
activity types and that they evolve independently. We 
search for the form of evolution that best reproduces 
the existing data. The new optimisation of the model 
parameters reproduces: (i) the number counts at 15, 24, 
60, 70, 90, 160, 170, and 850 /im, (ii) the known redshift 
distributions (mainly at 15 and 170 /im), (iii) the local 
luminosity functions at 60 and 850 /im, (iv) the CIB 
(from 100 to 1000 /mi) and its fluctuations (at 60, 100 
and 170 /im). 

Compared with the form of the model derived by 
Lagache et al. (2003), only a slight change of the co- 
moving luminosity density distribution is required (Fig. 
^), together with minor modifications to the starburst 
template spectra mainly between 12 and 30 /mi (Fig. |2J. 
The modified starburst spectra still reproduce the color 
diagrams of Soifer & Neugeubauer 1991 (12/25 versus 
60/100). We do not modify any other parameters in 
the model, nor the SED of the 'normal' galaxy template 
nor the SED of the starburst galaxy templates at longer 
wavelengths, nor the parametrisation of the local LF. 

PAHs radiate about 10% of the 1-1000 /im luminosity 
in a set of features concentrated mostly between 6 and 
15 /mi (Fig. I2J • This set of features is the strongest 

3 This assumption is based primarily on observations by ISO- 
CAM and ISOPHOT, but is confirmed by the first Spitzer studies 
of galaxy SEDs in the Lockman Hole and Groth Strip (Le Floc'h 
et al. 2004; Alonso Herrero et al. 2004) 
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TABLE 1 

1 cr c CONFUSION NOISE VALUES USING THE BEST CONFUSION 

ESTIMATOR OF DOLE ET AL. (2003), CONFUSION LIMIT Sn m AND 

THE VALUE OF q = Sii m /<7 C 



24 /im 



70 /tin 



160 /im 



1 <Tc 


8.0 /ijy 


^lim 


56 ,uJy 


q 


7.0 



0.47 mjy 
3.2 mJy 

6.8 



10.6 mjy 
39.8 mJy 

3.8 



a Using the source density criterion 

b In this case, the photometric and source density criteria agree 




Fig. 1. — Co-moving luminosity density distribution of the 
Lagache et al. (2003) model (dashed line) and the updated model 
presented here (continuous line). Also shown for comparison is the 
co-moving luminosity density distribution from all cases of Gispert 
et al. 2001 (crosses with error bars), together with the best fit 
passing through all cases (dot line). 



spectral concentration known in galaxy spectra. The 
features affect strongly the mid-IR number counts when 
the redshift brings them into the bandpass filter. For 
the 15 /im ISOCAM observations, the peak contribution 
arises when the 7.7 /Km maximum is close to the mid- 
dle of the band i.e. z-0.8. The ISOCAM 15 fjm counts 
showed the existence of a maximum in the redshift dis- 
tribution just at that redshift. Futhcrmore, it leads to 
a maximum in the counts (normalised to the Euclidian, 
e.g. on Fig. |3J, which can exist only when the IR energy 
output in the critical redshift range is dominated by a 
well defined luminosity range. Our model predicted a 
similar behaviour for the 24 /im band. 

3. RESULTS 

3.1. Source counts 

We show the comparison of the number counts at 15, 
60, 170 and 850 /um with the observations and as well as 
the new Spitzer 24 (im (Papovich et al. 2004), 70 and 
160 /itm (Dole et al. 2004a) counts in Fig. Eland Fig.H 
The updated model reproduces very well the counts at 15 
and 24 /im without any significant changes in the other 
number counts, as well as the CIB and its fluctuations 
and the redshift distributions of resolved sources at 15, 
60, 170 and 850 /im. The agreement between the ob- 
served counts and the model at 24 /xm is excellent. 

3.2. First interpretation of the Spitzer /MIPS results 




Lambda (/itn) 

Fig. 2. — Starburst model spectrum for L =3 10 11 Lq. Dashed 
line: Lagache et al. (2003), continuous line: the present updated 
version of the model. The spectrum has not been modified for 
A >30 /im. The biggest change is about a factor 2 (for a 3 10 11 
Lq galaxy) around 15 /tm. 
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Fig. 3. — Number counts at 15, 60, 170 and 850 /im (in Log) 
together with the model predictions (present work: continuous line, 
Lagache et al. 2003: dashed line). Data at 15 /tm are from Elbaz 
et al. (1999), at 170 /tm from Dole et al. (2001), at 60 /im from 
Hacking & Houck (1987), Gregorich et al. (1995), Bertin et al. 
(1997), Lonsdale et al. (1990), Saunders et al. (1990) and Rowan- 
Robinson et al. (1990) and at 850 /im from Smail et al. (1997), 
Hughes et al. (1998), Barger et al. (1999), Blain et al. (1999), 
Borys et al. (2002), Scott et al. (2002) and Webb et al. (2002). 



The consistency of the model from the mid-IR and 
far-IR to the submm indicates that the underlying as- 
sumptions are valid up to redshift 2.5. In particular, 
the very close agreement between the model and number 
counts at 15 and 24 /im implies that the PAH features 
remain prominent in the redshift range 0.5 to 2.5. Fur- 
thermore, the well-defined bump in the 24 /tm number 
counts (normalised to the Euclidian) at S~0.3 mJy (Fig. 
0J) implies at z~2 that Lg^m^ 3.2 10 11 L Q which corre- 
sponds to Li_iooo/im=3.5 10 12 L Q . These luminosities 
have to dominate the IR energy ouput at that redshift 
(as shown on Fig. UJ. As a comparison, an IR energy 
output dominated by L* galaxies at z~2 would maximise 
the 24 /im counts around 3 /iJy which is totally inconsis- 
tent with the observation since most of the CIB at 24 /im 



PAH contribution to the IR output energy at z~2 
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Fig. 4. — Number counts at 24, 70 and 160 £tm together with the model predictions (present work: continuous line, Lagache et al. 2003: 
dashed line). Data at 24, 70 and 160 fim are from Papovich et al. (2004) and Dole et al. (2004a). Notice that none of the observed source 
counts arc corrected for uncomplctcness. 



is resolved at flux densities brighter than ^60 /LtJy. 

Using the new model of IR galaxy evolution, Dole et 
al. (2004b) have made new determinations of the con- 
fusion limits for MIPS (Multi Band Imaging Photome- 
ter for Spitzer). A summary is provided in Table ^ 
Fig. [5] shows the updated redshift distributions for such 
confusion- limited surveys. At 24 /xm, the deepest sur- 
veys will probe star formation up to redshift 3. In the 
far-IR, MIPS surveys will probe the largely unexplored 
window l<z<2. 

In Fig El we show the different redshift contributions to 
the 24 /im number counts. The peak in the counts is 
dominated by galaxies with redshift between 1 and 2. 
There is also a non negligible contribution from z >2 
galaxies (30% at S=0.2 mjy). The integral of the source 
counts down to 60 /iJy gives 1.9 nW m~ 2 sr _1 , in ex- 
cellent agreement with Papovich et al. (2004). It repre- 
sents about 63% of the CIB at 24 /an . It should also be 
noted from Fig. [5] that the redshift distribution changes 
sharply for fluxes between 0.3 and 2 mjy. For sources 
weaker than 150 //Jy, the redshift distribution remains 
about constant. 

4. CONCLUSION 

We have updated the model of Lagache et al. (2003) to 
fit the new constraints from Spitzer data. We show that 
only a small change of the co-moving luminosity density 
distribution and of the starburst template spectra are 
required to fit all the observations from the mid-IR to 
the sub-millimeter. The agreement between the model 
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Fig. 5.— MIPS 24, 70 and 160 /im redshift distributions (in 
logarithmic scale). Solid lines: Shallow surveys (flux greater than 
0.3, 22.3 and 135.1 mjy at 24, 70 and 160 /im respectively), dashed 
line: Deep surveys (flux greater than 0.094, 6.8, 50.2 mjy at 24, 
70 and 160 (im respectively), Dotted line: Ultra Deep survey (flux 
greater than 0.06 mjy at 24 /an). 
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Fig. 6. — Redshift contribution to the number counts at 24 (im. 
The dot, dash, dash-dot, dash-3 dot, long-dash correspond to the 
number counts up to redshifts 0.3, 0.8, 1, 1.3 and 2 respectively. 
Data are from Papovich et al. (2004). 
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Fig. 7. — Co-moving evolution of the IR energy output per 
dex of luminosity (L dN/dLogL). The continuous, dashed, dashed- 
3dottcd lines are for z=0, 0.5 and 2 respectively. 



and all the data demonstrates that the integrated SEDs 



from galaxies still have prominent 6-9 /im emission 

from PAH features up to redshift 2.5. It clearly shows 

that the population of IR galaxies is not dominated by 

galaxies with featureless continuum spectra (such as 

AGN-type SEDs in the mid-IR). The fraction of the 

energy radiated by the PAHs is about one tenth of the 

total IR radiation locally and remains roughly constant 

between redshifts 0.5 and 2.5. 

Our initial model predicted that the IR output energy 

had to be dominated by galaxies with luminosities 

increasing rapidly with redshift (Li_iooopm— 3 10 12 L Q 

at z=2.5). This was needed mainly to account for the 

SCUBA counts (the total energy ouput at that redshift 

being mostly driven by the submillimeter extragalactic 

background level). The excellent agreement between 

the 24 [im counts and the model confirms strikingly 

the evolution of the IR luminosity function as being 

dominated in energy by ~3 10 11 L Q to ^3 10 12 L Q 

galaxies from redshift 0.5 to 2.5. 

We have for the first time with Spitzer an unbiased 

view of the IR Universe up to redshift 2.5. With all 

the mid-IR to sub-mm data, the model is now very well 

constrained up to z~2.5. Other predictions, together 

with the number counts of all the planned and future 

mid-IR to mm surveys, will be available on the web page: 

http : //www . ias . f r/PPERSD/glagache/act/gal$_$model . html. 
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